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Abstract. We carried out 1.25 pc resolution MHD simulations of 
the ISM - including the large scale galactic fountain ~ by fully tracking 
the time-dependent evolution of the magnetic field and the formation 
of shock compressed regions in a supernova-driven ISM. The simulations 
show that large scale gas streams emerge, driven by SN explosions, which 
are responsible for the formation and destruction of shocked compressed 
layers with lifetimes up to some 15 Myr. The T< 10^ K gas is dis- 
tributed in filaments which tend to show a preferred orientation due to 
the anisotropy of the flow induced by the galactic magnetic field. The 
simulations also show that the magnetic field has a high variability, it is 
largely uncorrelated with the density and it is driven by inertial motions. 
The latter is consistent with the fact that ram pressure dominates the 
flow for 10^ <T< 10^ K. For T> 10® K thermal pressure dominates, 
while for T< 10^ K (stable branch) magnetic pressure takes over. 
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1. INTRODUCTION 

Recent work by Kim et al. (2001) and Passot & Vazquez- 
Semadeni (2003) show somewhat contradicting results on the scaling 
of the magnetic field with density in the ISM. Kim et al. using a (200 
pc)'^ box centered in the Galactic midplane with periodic boundary 
conditions, driven by SNe at a rate of 12 times the Galactic value 
and using an uniform field strength of 5.8 /iG orientated along the 
x— direction, claim that the magnetic field scales as p'^-^ at densities 
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n > 1 cm^^, although their Fig. 2 shows an almost order of magnitude 
variation in the field for the same density. On the other hand Passot 
& Vazquez-Semadeni show that for small Alfvenic Mach number the 
magnetic field is uncorrelated with density, exhibiting a large scatter, 
which decreases towards higher densities. 

In this paper we investigate the variability of the magnetic field 
in the Galactic disk and its correlation with the density in global 
MHD simulations of the ISM, driven by SNe occuring at the Galactic 
rate. Other important issues like the volume filling factors of the ISM 
"phases", the dynamics of the galactic fountain, the conditions for 
dynamical equilibrium and the importance of convergence of these 
results with increasing grid resolution have been treated elsewhere 
(Avillez 2000, Avillez & Breitschwerdt 2004a, b). 

2. MODEL AND SIMULATIONS 

We ran kpc-scale MHD simulations of the SN-driven ISM on a 
cartesian grid of < {x,y) < 1 kpc size in the Galactic plane and 
— 10 < < 10 kpc into the halo with a finest adaptive mesh re- 
finement resolution of 1.25 pc, using a modified version of the 3D 
model of Avillez (2000). The new model allows for the motion of 
OB associations and includes magnetic fields (Avillez & Breitschw- 
erdt 2004b). At the beginning of the simulations the uniform field 
components along the three axes are given by (i?„,o("^(^)/"^o)^^^, 0, 0), 
where Bufi = 3 is the field strength, n{z) is the number den- 
sity of the g function of distance from the Galactic midplane 
and no = lcm~^ is the average midplane density. The random field 
component is set to zero in the beginning of the simulations. This 
component is built up during the first millions of years of evolution 
as a result of turbulent motions, mainly induced by SN explosions. 

3. RESULTS 

3.1. Global Evolution 

The initial evolution of the magnetized disk is similar to that 
seen in purely HD runs (Avillez 2000, Avillez & Breitschwerdt 2004a), 
that is, the initially stratified distribution does not hold for long as 
a result of the lack of equilibrium between gravity and (thermal, ki- 
netic and turbulent) pressure during the "switch-on phase" of SN 
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Left panel: mavillezJigla, Right panel: mavillezJiglb 



Fig. 1. Density and magnetic field distribution in the Galactic midplane 
after 374 Myr of evolution. The resolution of the finest AMR level is 1.25 
pc. 

activity. As a consequence the gas in the upper and lower parts of 
the grid collapses onto the midplane, leaving low density material 
in its place. However, in the MHD run it takes a longer time for 
the collapse to be completed as a result of the opposing magnetic 
pressure and tension forces. As soon as the system has collapsed 
and enough supernovae have gone off in the disk building up the 
required pressure support, transport into the halo is not prevented, 
although the escape of the gas takes a few tens of Myr to occur. The 
crucial point is that a huge thermal overpressure due to combined 
SN explosions can sweep the magnetic field into dense filaments and 
punch holes into the extended warm and ionized Hi layers. Once 
such pressure release valves have been set up, there is no way from 
keeping the hot over-pressured plasma to follow the pressure gradient 
into the halo. As a consequence the duty disk-halo-disk cycle of the 
hot gas is fully established, in which the competition of energy input 
and losses into the ISM by SNe, diffuse heating, radiative cooling 
and magnetic pressure leads the system to evolve into a dynamical 
equilibrium state within a few hundred Myr. This time scale is con- 
siderably longer than that quoted in other papers (e.g., Korpi et al. 
1999, Kim et al. 2001), because in these the galactic fountain has 
not been taken into account. 

Fig. 1 shows slices of the 3D data cube of the density and mag- 
netic field distributions in the Galactic midplane. The highest den- 
sity gas tends to be confined to shocked compressed layers that form 
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in regions where several large scale streams of convergent flow (driven 
by SNe) occur. The compressed regions, which have on average life- 
times of 10-15 Myr, are filamentary in structure, tend to be aligned 
with the local field and are associated with the highest field strengths. 
The formation time of these high density structures depends on how 
much mass is carried by the convergent flows, how strong the com- 
pression is and on the rate of cooling of the regions under pressure. 



3.2. Field Variability and Dependence with Density 
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During the evolution of the 
system, thermal and dynami- 
cal processes broaden the dis- 
tribution of the field strength in 
such a way that after the global 
dynamical equilibrium has been 
set up the field strength in the 
disk spans two orders of magni- 
tude from 10-^ to 10-^ G (Fig. 
2). The figure shows a large 
scatter in the field for the same 
density suggesting that the field 
is being driven by the inertial 
motions, rather than it being 
the agent determining the mo- 
tions. In the latter case the 
field would not be strongly dis- 
torted, and it would direct the 
motions predominantly along 
the field lines. The high field 
variability is also seen in the 

right panel of Fig. 1, which shows a highly turbulent field, that seems 
to be uncorrelated with the density, and thus, the classical scaling 
law B ~ p", with a = 1/2, according to the Chandrasekhar-Fermi 
(CF) model (1953) will not hold. It should be kept in mind that in 
CF it was assumed that the field is distorted by turbulent motions 
that were subalfvenic, whereas in our simulations in addition both 
supersonic and superalfvenic motions can occur, leading to strong 
MHD shocks. It now depends if the density fluctuations are primar- 
ily caused by thermal and/or ram pressure fluctuations and there- 
fore being uncorrelated with the magnetic fleld, or if MHD waves 



Fig. 2. Scatter plot of B vs /? for the 

T< 10''^ (black), 10^ <T< 10^ (green), 
10^ <T< lO^-^ K (blue) and T> lO^-^ 
K (red) regimes at 400 Myr of disk evo- 
lution. 
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Fig. 3. The left panel shows the average magnetic (red), thermal (black) 
and ram (green) pressures as functions of temperature at 400 Myr. The 
right panel shows the variation of with Pt/^/T for T< 10^ K gas in the 
disk. The 2a parameters are the slopes of the straight lines and correspond 
to the relation P_b cx (Pj/^/T)'^" if a i? oc law is assumed. 



are the driving agent and provide a coupling between matter and 
field. Therefore, in general < a < 1 would be expected. However, 
it should be noted that in reality heating and cooling processes, and 
even magnetic reconnection could induce further changes. Observa- 
tionally there seems to be some evidence from measuring magnetic 
field strengths in the cold neutral medium, that B and p can be 
largely uncorrelated (Troland & Heiles 2001). 

The left panel of Fig. 3 shows that the T< 10^ K gas has Pb > 
Pram > Pth, demonstrating that magnetically dominated regions do 
exist, while the T> 10® K gas has the highest thermal pressure and 
the lowest magnetic pressure. For T< 10^ K the relation B oc p°' does 
not hold either, as a varies between -0.006 and 0.085 (right panel 
of Fig. 3) suggesting that the thermal and magnetic pressures are 
independent. This result is largely consistent with the almost zero 
variation of magnetic pressure with temperature seen in the left panel 
of Fig. 3. For 10^ <T< 10® K ram pressure determines the dynamics 
of the flow, and therefore, the magnetic pressure does not act as a 
restoring force (Passot & Vazquez-Semadeni 2003) as it was already 
suggested by the lack of correlation between the field strength and 
the density. Fig. 3 also shows that the basic assumption of energy 
equipartition, made in the paper by CF (1953) in order to calculate 
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the magnetic field in the spiral arm, is clearly not fulfilled. 



4. DISCUSSION 

The dynamical picture that emerges from these simulations is 
that thermal pressure gradients dominate mostly in the neighbor- 
hood of SNe, which drive motions whose ram pressures are dominant 
over the mean thermal pressure (away from the energy sources) and 
the magnetic pressure. The magnetic field is only dynamically im- 
portant at low temperatures, but can also weaken gas compression 
in MHD shocks and hence lower the energy dissipation rate. The 
thermal pressure of the freshly shock heated gas exceeds the mag- 
netic pressure by usually more than an order of magnitude and the 
B-field can therefore not prevent the fiow from rising perpendicular 
to the galactic plane. Thus hot gas is fed into the galactic fountain 
at almost a similar rate than without field. 

The present simulations do not include self-gravity, which could 
have an important effect on high density regions, although on the 
verge of Jeans instability these would decouple from the ambient ISM 
fiow. Moreover heat conduction has been neglected in the view of it 
being second order in comparison to the dominance of bulk motions 
and turbulent mixing. In the spirit of building a bottom-up model of 
the ISM further components (e.g., cosmic rays) and processes (e.g., 
non-equilibrium cooling) will be included in future simulations. 
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